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Abstract. For models which exhibit a continuous phase transition in the 
thermodynamic limit a numerical study of small systems reveals a non- 
monotonic behaviour of the microcanonical specific heat as a function of the 
system size. This is in contrast to a treatment in the canonical ensemble where 
the maximum of the specific heat increases monotonically with the size of the 
system. A phenomenological theory is developed which permits to describe this 
peculiar behaviour of the microcanonical specific heat and allows in principle 
the determination of microcanonical critical exponents. 

PACS numbers: 05.50.-|-q,64.60.-i,65.40.Gr 
1. Introduction 

In recent years numerous studies investigated the possible differences between 
the microcanonical and the canonical treatment of a given system. It is now well 
accepted that in various cases the microcanonical and the canonical ensemble are 
not equivalent [Ill21EllllElinilZIIHl- For short-range interactions, the equivalence 
of the two ensembles holds in the infinite volume limit, but this is not the case 
in finite systems. For long-range interactions, as encountered for example in 
gravitational systems, the two ensembles remain inequivalent even for infinite 
systems. 

Clearly, this inequivalence in finite systems makes the microcanonical 
analysis of possible signatures of phase transitions an important issue j21 El UHl HI 
ElEj- For discontinuous phase transitions, the microcanonical analysis reveals 
typical small system signatures, as e. g. a back-bending of the caloric curve or 
the appearance of a negative specific heat. Negative heat capacities have indeed 
been measured in recent experiments on nuclear fragmentation ^21 and on the 
melting of atomic clusters ^1]. Similarly, intriguing features are also revealed in 
the microcanonical analysis of small systems which exhibit a continuous phase 
transition in the thermodynamic limit. Indeed, typical features of symmetry 
breaking, as e. g. the abrupt onset of a non-zero order parameter when the 
(pseudo-) critical point is approached from above or a diverging susceptibility. 
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turn up already for finite systems |9j . This is in contrast to the canonical ensemble 
where singularities appear exclusively in the thermodynamic limit. 

The fundamental quantity in a microcanonical analysis is the density of 
states or, equivalently, the microcanonical entropy. All relevant quantities can 
indeed be expressed by partial derivatives of the microcanonical entropy. For 
example, the susceptibility is proportional to the inverse of the curvature of the 
entropy surface. It is the existence of a point with vanishing curvature that is 
responsible for the divergent susceptibility observed in finite systems which have 
a continuous phase transition in the infinite volume limit. 

In the present work we examine more closely the finite-size behaviour of 
the microcanonical specific heat in different classical spin systems. For systems 
with a continuous phase transition one expects that the maximum of the specific 
heat increases with increasing system size. This is indeed observed in the 
microcanonical analysis for not too small systems. For small system sizes, 
however, we observe a non-monotonic behaviour as the maximum of the specific 
head first decreases for increasing system sizes. This is again a property of the 
entropy surface as the microcanonical specific heat can be exclusively expressed 
by energy derivatives of the microcanonical entropy. In order to account for this 
peculiar behaviour we develop a phenomenological theory based on the analyticity 
of the entropy surface of finite systems. 

The paper is organized as follows. In Section 2 we first discuss the definition 
of the temperature in the microcanonical ensemble. In the microcanonical 
ensemble various definitions of the temperature are possible, the different 
expressions becoming equivalent in the thermodynamic limit. The microcanonical 
specific heat, based on the expressions for the temperature, is the subject of 
Section 3. Numerical results obtained for two- and three-dimensional Ising 
models as well as for the two-dimensional three-state Potts model reveal a non- 
monotonic behaviour of the specific heat for increasing system sizes. The finite- 
size behaviour of the specific heat of microcanoncial systems is considered from 
a phenomenological point of view in Section 4 where a finite-size scaling theory 
is developed which explains the peculiar behaviour of the microcanonical specific 
heat of small systems. Finally, Section 5 gives our conclusions. 

2. Temperature in the microcanonical ensemble 

The density of states is the starting point for the statistical description of 
thermostatic properties in the different ensembles. For a magnetic system that is 
isolated from any environment the proper natural variables are the energy E and 
the magnetisation M. The corresponding characteristic function of the isolated 
system is the microcanonical entropy 

S{E,M,L-^) =\nn{E,M,L-^) (1) 
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where Q denotes the degeneracy of the macrostate {E, M) and L is the hnear 
extention of the system. Here and in the following units with ks = 1 are 
used. The microcanonical analysis of finite classical spin systems starts from 
the microcanonical entropy density 

s(e, m, L-') := ■^S{L'^e, L'^m, L-^) (2) 

of a system in d dimensions with N = L'^ spins, where e = E/N denotes the 
energy density and m = M/N the magnetisation density. In the following the 
dependence on the system size is suppressed in order to improve readability. 

Before investigating the microcanonical specific heat, we first have to 
discuss the definition of the temperature for finite systems in the microcanonical 
ensemble. In the thermodynamic limit canonically defined physical quantities and 
the corresponding microcanonical quantities have to become identical for systems 
with suitably short range forces. However, this requirement does not yield an 
unambiguous definition of the microcanonical temperature, leading to different 
physically plausible definitions in finite systems which all become equivalent in 
the thermodynamic limit. 

The starting point is the canonical partition function 

Z{/3,h) = L'^'^ J de ^ dm exp |L^(s(e, m) - /3e + /3/im)| (3) 

which is the Laplace transform of the density of states. Here the inverse canonical 
temperature (3 and the applied magnetic field h are external parameters which 
are imposed on the system by its environment. The canonical temperature 
and external field are denoted by a tilde in order to avoid any confusion 
with the microcanonical temperature and field defined in the following. The 
integral (jHl) can be evaluated in the limit L ^ oo by means of the Laplace 
method. For a given inverse temperature jS and external magnetic field h the 
dominant contributions to the integral arise from the maximum of the argument 
g{e,m) = s{e,m) — (3e + (3hm. The equations deg = and dmg = suggest the 
following definitions of the inverse microcanonical temperature: 

/3(e,m) = des{e,m) (4) 

and of the microcanonical magnetic field: 

P{e,m)h{e,m) = —dms{e,m). (5) 

Here and in the following the notation dx is used for the partial derivative d/dx. 
The inverse microcanonical temperature (3 and magnetic field h are conjugate 
variables of the natural variables e and m of the microcanonical approach 
and consequently depend on these. The definition of the microcanonical 
temperature surface leads to the foUwing definition of the temperature in 
equilibrium. Consider the spontaneous magnetisation msp(e) of the magnetic 
system for a given energy e that is defined by the condition h{e, m) = The 
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temperature of the magnetic system in equilibrium is then obtained by evaluating 
the inverse temperature surface /5(e, m) at the equilibrium macrostate (e, msp(e)): 

(3^(e) := /3(e,msp(e)) = des{e,m)\m=msp{e)- (6) 

This definition of the inverse temperature ensures ensemble equivalence between 
the canonical and microcanonical description, as can be seen using the Laplace 
method in the asymptotic limit L oo. For finite L, however, the 
exponential in ^ cannot be approximated by the quadratic term of the Taylor 
expansion only. Higher order terms are necessary which render the integrand 
asymmetric. Consequently, the canonical mean values are shifted from the 
associated maximum of the entropy surface leading to the inequivalence of the 
canonical and the microcanoncial ensemble for finite system sizes. 

We pause here for a moment to recall that the spontaneous magnetisation 
^sp(e) of a finite microcanonical magnetic system exhibits features which are 
typical of phase transitions. The spontaneous magnetisation of the Ising model 
in dimensions d > 2, for example, is zero above a well-defined transition energy 
Cpc and becomes non-zero below Cpc. Close to this pseudo-critical energy the 
variation of the spontaneous magnetisation as a function of the deviation of e 
from Cpc is described by a square root function [HITTj. This classical behaviour 
has its origin in the analyticity of the entropy surface for all finite systems 
The appearance of a non-zero spontaneous magnetisation reflects the spontaneous 
breakdown of the global symmetry of the system and may be regarded as a 
precursor of the critical point of the infinite system jHl HH HZI. Note that the 
specific entropy Soo(e, m) in the thermodynamic limit is a concave function of its 
variables. In finite systems, however, this is not compulsory so that two maxima 
of the entropy can appear at non-zero magnetisations for a given energy. 

Coming back to the canonical ensemble we remark that in absence of a 
magnetic field the partition function Q simplifies to 

Z{(3,h = 0) = L"^ j de (^L'^ j dm exp{L°'s(e,m)}^ exp{-L^/3e},(7) 

which leads to the definition 

exp{L'^SR(e)} = L'^ j dm exp {L^s(e,m)} (8) 

of the reduced (specific) entropy SR(e). In the limit of large system sizes the 
dominant contributions to the integral ((Zj) arise from the energy defined by the 
maximum of the argument ^^(e) — f3e for a given inverse canonical temperature 
/3. This suggests the following alternative definition of an inverse (reduced) 
microcanonical temperature, namely 

/5R(e) = ^SR(e). (9) 

The thermal properties of the system are now obtained from the entropy function 
SR(e) rather than from the full entropy surface s{e,m) depending on both the 
energy and the magnetisation. 
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To conclude this section the interrelation between the inverse temperatures 
I3e{g) and Pji{e) is briefly considered. In the asymptotic limit L ^ oo the integral 
(jHl) is dominated by the entropy s(e,msp(e)), evaluated at the spontaneous 
magnetization msp(e), as can again be seen by using the Laplace method. 
Therefore, the entropy s^{e) is given by s(e,msp(e)) for asymptotically large 
system sizes L and one gets 

/3R(e)'^~°°^3(e,m,p(e)). (10) 

Carrying out this differentiation we obtain the relation 

L d d d 

Pnie) — s(e,m)|^=^,p(e) + — s(e, m)|^=^,p(e)— msp(e). (11) 

As dms{e,m) is zero at (e,msp(e)) the second term vanishes and one is left with 
Pr{^) ~ Pe{g)- The full entropy surface s(e, m) and the reduced entropy function 
spi,(e) will therefore lead to the same equilibrium temperature in the asymptotic 
limit L — s> oo. For finite L, however, /?E(e) is significantly different from 



3. Microcanonical specific heat 
3.1. General discussion 

Once the inverse temperature /5 of a microcanonical system is evaluated — here /3 
may be /3e or /5r — one can calculate the specific heat which is generally defined 
by c = du/dT with u and T = 1/(3 being the energy and the temperature of the 
system. For the microcanonical specific heat as a function of the energy of the 
system this gives 

The discussed ambiguity in the definition of the microcanonical temperature leads 
also to different expressions for the specific heat in finite systems. However, they 
converge towards the same limit function in the thermodynamic limit. 

In the following we discuss the finite-size behaviour of the specific heat arising 
from the temperature (3ji{e) in different classical spin models (from now on we 
drop the subscript R in order to avoid unnecessary notation). Note that this is 
the definition of the specific heat that is the most relevant for experiments where 
usually the energy is considered as the unique natural variable corresponding 
to systems to which no external field is applied. Specifically, we study three 
models undergoing a continuous phase transition in the thermodynamic limt: 
the two- and and the three-dimensional Ising model as well as the three-state 
Potts model in two dimensions. The nearest neighbour Ising model is defined by 
the Hamiltonian 

7i = - ^cTiO-j, (13) 
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where the summation over nearest neighbour pairs is indicated by {i,j) and the 
spin (jj at site i can be in the states (Tj = ±1. In the present study Ising models 
defined on the square and on the cubic lattices are considered. The three-state 
Potts model is a generalisation of the Ising model where the Potts spins cxj take 
on the values 1,2,3. The Hamiltonian is given by 

n = -Y,6,^,^^, (14) 

where = 1 when the spins located at the neighbouring sites i and j have the 
same value and zero otherwise. For the Potts model we only consider the square 
lattice. 

In finite systems the appearance of a continuous phase transition in the 
thermodynamic limit is signaled by a maximum in the specific heat which becomes 
more and more pronounced when the system size is increased. This behaviour 
of the specific heat is due to a maximum of the second derivative of s(e) which 
is negative everywhere and tends to zero for increasing system sizes from below. 
The position of the maximum of the second derivative of s(e) defines a pseudo- 
critical energy Cpc of the finite system. At the same time the microcanonical 
inverse temperature /3pc := /?(epc) evaluated at the energy Cpc converges towards 
the critical value 1/Tc when L tends to infinity. 



des{e) dees{e) 




Figure 1. Schematic depiction of the first two derivatives of the reduced 
entropy near the pseudo-critical energy of a finite system. The pseudo- 
critical energy Cpc corresponds to an inverse temperature /3pc (left). The 
derivative deeS has a maximum at Cpc (right). 

The behaviour just described, which is schematically sketched in figure ^ 
is indeed observed in the different models for not too small system sizes. For 
very small systems, however, our numerical results reveal an unexpected non- 
monotonic behaviour of the specific heat, as discussed in the next subsection. 
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Figure 2. Evolution of the coefficient ^2(0 finite three- 

dimensional Ising systems with periodic boundary conditions (a) and 
open boundaries (b). The modulus of the coefficient ^2(0 of the system 
with periodic boundary conditions has a plateau and becomes smaller 
again for small systems. In the thermodynamic limit I — > both curves 
have to extrapolate to A2 = 0. 



3.2. Numerical results 

In this subsection the specific heat of finite Ising and Potts systems is investigated 
numerically for both periodic and open boundary conditions. To obtain the 
numerical data we used a recently proposed very efficient method for the direct 
computation of the density of states JT] . Specifically, we discuss in the following 
the value deeS^Cpc, L~^) of the second derivative of the entropy evaluated at the 
pseudo-critical energy Cpc- For later convenience this value is denoted by ^2(0 
where / := is the inverse system size: 

^2(0 = '9ees(epc,0- (15) 

The coefficient ^2(0 of finite three-dimensional Ising systems is shown in 
figure 121 With periodic boundary conditions the coefficient A2{1) shows a back- 
bending as its modulus first increases for increasing / (i. e. decreasing system 
sizes) and then decreases for very small systems, see figure This intriguing 
and unexpected back-bending is not observed in the system with open boundaries. 

Similarly, the coefficient A2{1) of the two-dimensional Ising model with 
periodic boundary conditions also exhibits this back-bending for very small 
system sizes, whereas again no back-bending is observed for open boundaries, 
see figure El In case of the systems with linear extensions L = 8, 10 and 32 and 
periodic boundaries the numericaly determined data can be compared to exactly 
computed data [IHl CHI 120] ■ This is also indicated in figure El 

Naturally, the back-bending of the coefficient ^2(0 directly affects the 
behaviour of the microcanonical specific heat of small systems as can be seen from 
equations (fT3|) and (fT^ . Indeed, the maximum of the specific heat first decreases 
with growing system size before increasing again, thus yielding a divergence in the 
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Figure 3. Evolution of the coefficient ^2(0 finite Ising models in two 
dimensions with periodic boundary conditions (a) and free boundaries 
(b). The open circles in (a) display exactly evaluated data. The 
coefficient ^2(0 bends back for small systems only for periodic boundary 
conditions. For the system with open boundaries a logarithmic plot 
shows that the coefficient ^2(0 evolves logarithmically for large systems 
(see also section 4). 



thermodynamic limit. This decrease of the specific heat of small microcanonical 
systems is displayed in figure|3]for the three-dimensional Ising model with periodic 
boundary conditions. It is worth noting that such a peculiar behaviour of the 
specific heat of small systems is not observed in the canonical ensemble (see, 
e.g. [mini 122) ■ 




Figure 4. Evolution of the microcanonical specific heat for small three- 
dimensional Ising systems with periodic boundary conditions. The 
maximum decreases for increasing system size. 



Finally, figure El displays the evolution of the coefficient ^2(0 fo'^ finite 
two-dimensional three-state Potts models for both periodic and open boundary 
conditions. The back-bending is strongly pronounced for periodic boundaries 
and is in this case also visible, but less developed, for free boundaries. For large 
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systems the coefficient ^2(0 eventually approaches zero reflecting the appearance 
of a continuous transition in the infinite system. 
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Figure 5. Evolution of the coefficient ^2(0 finite three-state Potts 
models in two dimensions with periodic boundary conditions (a) and 
free boundaries (b). For both boundary conditions the coefficient ^2(0 
bends back for small systems. 



4. Phenomenological theory for finite systems 

In this section the behaviour of the specific heat of finite microcanonical systems is 
investigated from a phenomenological point of view. As discussed in the following 
this leads to a theoretical description that accounts for the peculiar behaviour of 
the microcanonical specific heat described in the previous section. 

The specific heat of the infinite Ising or Potts systems diverges at the critical 
point ec,oo- For a system with a power law singularity the specific heat has the 
form 

Coo(e) ~ |e - Ccool""^ (16) 

in the vicinity of the critical point, where denotes the microcanonical critical 
exponent. The specific entropy of the infinite system contains a singular part 
that is a generalised homogeneous function in the vicinity of ec,oo characterised 
by the degree of homogeneity . The microcanonical critical exponent is 
related to by 

= = (17) 

I — a 

with a being the critical exponent of the canonical specific heat ^9;. Similarly, 
the microcanonical critical exponent i'^ of the correlation length can be expressed 
as 

"'-^-T^ (18) 

where the dimensionality of the system is again denoted by d. 
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The discussion of finite-size scaling relations of the specific heat starts from 
the decomposition 

s{e,l) = Sr{e,l) + s,{eJ) (19) 

of the entropy of a finite system into a regular and a singular part. The deviation 
of the energy from the pseudo-critical energy is denoted by e := e — Cpc, where / 
again is the inverse system size. The singular and regular parts of the entropy of 
the finite systems are chosen to approach the corresponding singular and regular 
part of the entropy of the infinite lattice: 

limss/r(e,0 = Soo,s/r{e)- (20) 

Note that the singular part of the entropy of a finite system is an analytic function 
due to the analyticity of thermodynamic potentials of finite systems. The singular 
part Ss{e,l) is assumed to obey the scaling assumption 

s,iEj) = jS,i\''^e,\'^''l) (21) 

with a positive re-scaling factor A and the degree of homogeneity discussed 
above. The ansatz ()2ip for the finite-size behaviour of the singular part of the 
microcanonical entropy of finite systems does not account for additional finite- 
size corrections that arise from the contributions of irrelevant scaling fields. The 
quahtative picture that is developed in the following can be extended to include 
those contributions as well. Differentiating the finite-size scaling assumption ^T^i 
with respect to the re-scaling factor A and setting A = 1 afterwards gives rise to 
the differential equation 

Ss{e, I) = a^ed^s^{e, + ^ idis^^e, I) (22) 

for the singular part of the reduced entropy. 

The entropy of a finite system is analytic |TH] and therefore Sg can be 
expanded with respect to the energy deviation e, yielding the series expansion 

oo _ 

s,{e,l)=Bo{l)+y2-B^{l)e\ (23) 

n=l 

From equation (j22j) the differential equation 

Bn{l) = na,Bn{l) + ^ ldiBn{l) (24) 

is obtained for the expansion coefficients Bn{l), n = 0,1,... This differential 
equation has the solution 

Bnil) = (25) 

where the Bn'' are size-independent coefficients. Similarly, the regular part of the 
reduced entropy can be expanded into the series 

oo ^ 

s,(£,0=Co(/) + V-C„(/)£". (26) 
^ — ^ n 

n=l 
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As the regular part of the entropy of the infinite system is also analytic, it is 
natural to assume that the coefficients are regular functions in I so that the 
Cn(0) are the expansion coefficients of the regular part of the entropy of the 
infinite system. Mathematically speaking, the limiting procedure / — > and the 
summation in ()26|) can be interchanged. This assumption is not possible for the 
singular part whose limit in the infinite system is non-analytic. 

Taking everything together we obtain that the expansion of the total entropy 
of a finite system is of the form 

oo ^ 

s{e, I) = Soil) + [B^^^l'^'-'^''^^ + Cn{l)) (27) 

n=l 

To proceed further let us consider the vicinity of the pseudo-critical point Cpc of 
the finite system of inverse length I. From a Taylor expansion we obtain (with 

e = e- Cpc) 

s(£, /) = So(/) + /?pc£ + \A2{l)e^ + \AS)e' + • • • (28) 
where the coefficients of the second and fourth order terms, 

^2(0 = 5^/- +C2(0 (29) 

and 

Ai{l) = Bfl'^+CS). (30) 

involve the critical exponents of the microcanonical system. 

The coefficient ^2(0 of second degree term is of particular interest as it 
describes the evolution of the microcanonical specific heat at the pseudo-critical 
energy e^c- Indeed, the curvature at Cpc as a function of the system size is given 
by 

d,,s{e = 0, /) = A2{1) = Bfh^ + C2{1) < 0. (31) 

As the function C2{1) is regular and has to vanish in the thermodynamic limit 
(Z — > 0) in order to produce a diverging specific heat at the critical point ec,oo of 
the infinite system, it has to be of the form 

C2(/) = Vil + ]^V2f + . . . (32) 

for small /. For a continuous phase transition the coefficient ^2(0 is negative 
for all inverse system sizes / (see equations and (fT^ ). therefore, the 

coefficient B^^ is also negative as it is the dominating one for the asymptotic 
limit of vanishing /. However, the sign of the coefficient C2{1) is not further 
restricted. The possible evolutions of the coefficient ^2(0 depending on the sign 
and variation of C2{1) are schematically shown in figure IHl A back-bending of 
the function A2{1) for decreasing system sizes can be caused by a large enough 
positive coefficient C2(/). The resulting minimum in A2{\) has the consequence 
that the maximum of the specific heat of small systems decreases with increasing 
systems size L and increases again in the limit of large systems. This is exactly 
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what we observe numerically. Thus, the phenomenological viewpoint developed 
in this section accounts for the peculiar behaviour of the specific heat of small 
microcanonical systems reported in section EiSl 



A2{1) A2{1) 




Figure 6. Schematic discussion of the possible size-dependence of the 
expansion coefficient ^2(0 (solid lines) of the reduced entropy function. 
Shown are two cases, namely a negative (left) and positive (right) 
contribution C2(/). For small / the variation is determined by the sum 
of the singular contribution B2{1) governed by the ratio cte/vs (dashed 
lines) and the contributions from C2{1) (dotted lines). For a rapidly 
growing positive contribution C2{1) the coefficient can have a minimum 
(right). 

Finally, let us note that the consideration of the limit / — >■ (i. e. L — > oo) 
allows in principle the determination of the ratio ae/z/^ (see for another 
recent discussion of this point). In the limit of large systems the evolution of the 
coefficient ^2(0 a function of the inverse system size / is governed by the ratio 
OL^jv^. In figure [7| a double-logarithmic plot of the coefficient |v42(/)| is shown for 
the Potts system with periodic boundary conditions. The data for large systems 
is in good approximation described by a straight line. The slope of this line is 
an estimate of the exponent ratio a^jv^. From our data we obtain the slope 
0.42 ± 0.03 which has to be compared with the exactly known value 2/5. The 
evolution of the coefficient ^2(0 large system sizes is indeed determined by the 
critical exponent ratio a-i^jv^. Hence, the microcanonical analysis of the evolution 
of physical quantities of finite systems allows, in principle, the determination of 
the true critical exponents characterising the critical behaviour of the infinite 
system (see ^] for a discussion of how to determine the order parameter 
critical exponent directly from the density of states of small systems). 

The picture is somehow different for the systems with open boundary 
conditions. There, the asymptotic regime governed by the exponent a^jv^ is 
not yet reached for the considered systems sizes [L < 84). This probably has its 
origin in the finite-size contributions of the free surfaces [21] which strongly affect 
the behaviour of small systems with open boundaries. In order to determine the 



Microcanonical specific heat 



13 



d -2.1 




-3 -2.5 



HI) 

Figure 7. Double- logarithmic plot of the coefficient 1^2(01 of the three- 
state Potts model in two dimensions with periodic boundary conditions. 
For large systems the data is approximately given by a straight line. The 
linear regression for the data of the largest systems {L between 48 and 
84) gives the slope 0.42 ± 0.03 (dashed hue). 



ratio ae/i'e in this case much larger system sizes have to be investigated. 

The phenomenological picture developed so far applies only to systems with 
an algebraically diverging specific heat in the thermodynamic limit. For systems 
with a logarithmically diverging specific heat (as encountered for example in 
the two-dimensional Ising model), the canonical and microcanonical critical 
exponents are identical. This suggests that the finite-size behaviour of the 
various physical quantities at the pseudo-critical point should be described by 
the same asymptotic law in terms of the system size in both ensembles. As 
the microcanonical specific heat at Cpc is basically given by the inverse of the 
coefficient A2{1), the asymptotic behaviour of the canonical specific heat 
suggests the form 

1M2(0 = In L + 6*2(0 . (33) 

Plotting 1/^2 (/) against InL as done in Figure for the two-dimensional Ising 
model with open boundaries shows that the asymptotic law (p^ indeed holds for 
large system sizes. 

To conclude the phenomenological considerations of this subsection a short 
remark about the corrections to scaling due to irrelevant scaling fields must be 
added. These additional correction terms are non-analytic in the inverse system 
size / and alter therefore the size- dependence of the coefficient B2{1). Denoting 
the non-integer exponent of the correction to scaling term by the expression 
of ^2(0 fo'^ ^ system with an algebraically diverging specific heat is given by 
(compare relation (jSU) 

A2{1) = 5f (1 + 61/"^ + ...) + C2{1). (34) 

Here C2{1) is again the correction term that arises from the regular part of the 
entropy. Note that a negative coefficient hi with suitably large modulus may also 
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cause the possible back-bending of the coefficient ^2(0 small system sizes. 
5. Conclusions 

Precursor effects of phase transitions can be very different in the microcanonical 
and in the canonical treatment of finite systems. The best known example 
is the appearance of a negative microcanonical specific heat in small systems 
that announces a discontinuous phase transition. But typical features are also 
encountered in the microcanonical ensemble in cases where a continuous phase 
transition takes place in the thermodynamic limit, the most intriguing being a 
divergent susceptibility already present in finite systems. 

In this work we have shown that the microcanonical specific heat also shows 
a peculiar behaviour for small systems that undergo a continuous phase transition 
in the thermodynamic limit. The observed initial decrease of the specific heat 
for increasing system sizes has to be compared to the behaviour in the canonical 
ensemble where a monotonic increase of the maximum of the specific heat is 
encountered. 

We have presented a phenomenological finite-size scaling theory that permits 
to explain this peculiar behaviour. This theory, which is based on the analyticity 
of the microcanonical entropy surface, uses as a variable the distance to the 
pseudo-critical point Cpc of a given finite system. This unusual ansatz has allowed 
us recently to extract the order parameter critical exponent directly from the 
density of states of small systems ^2j. The phenomenological finite-size scaling 
theory should therefore be viewed in the broader context of deriving a finite-size 
scaling theory in the microcanonical ensemble. 

There do exist some earlier attempts at a microcanonical finite-size scaling 
theory. A finite-size scaling theory for a microcanonical ensemble with the energy 
as its only natural variable was also formulated in (SHI- However, that work 
is based on a definition of the microcanonical entropy of finite systems that is 
different from the definition used in our work. In fact, the definition of the 
entropy used in [2E] has a major disadvantage. It is well known that the various 
statistical ensembles can be formulated in a unified way in terms of the extremal 
properties of Boltzmann's eta-function. These extremal properties have to be 
worked out under certain subsidiary conditions which are related to the way how 
the system is coupled to its environment in the different ensembles. This unified 
point of view is, however, not possible for the microcanonical ensemble considered 
in 

Microcanonical finite-size scaling relations were also considered in [HI EZ! 
for the whole entropy surface s{e,m, L~^). In those works the analysis of the 
entropy surface s{e,m, L~^) was carried out with respect to the transition point 
of the infinite system. This is different in the considerations of the present work 
where the relative deviation from the finite system transition point has been 
investigated. Microcanonical finite-size scaling relations were also investigated 
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in j2Hl- In that work the microcanonical quantities were basically defined as 
expectation values with respect to the microcanonical probability distribution 
Pe{M) ~ Q{E,M). The microcanonical quantities analysed in the present work 
are defined in a conceptionally different manner. 

Finally, let us note that in experiments on nuclear systems or atomic clusters 
knowledge of the infinite system is usually not available. Therefore, our scaling 
theory involving only quantities of the finite system considered seems to be the 
most appropriate for describing this kind of experiments. 
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